Numerical simulation of dark lanes in post— flare supra— arcade 
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We integrate the MHD ideal equations to simulate dark void sunwardly moving structures in 
post-flare supra-arcades. We study the onset and evolution of the internal plasma instability to 
compare with observations and to gain insight into physical processes and characteristic parameters 
of these phenomena. The numerical approach uses a finite-volume Harten-Yee TVD scheme to 
integrate the 1D^ MHD equations specially designed to capture supersonic flow discontinuities. 
The integration is performed in both directions, the sunward radial one and the transverse to the 
magnetic field. For the first time, we numerically reproduce observational dark voids described 
in Verwichte et al. (2005). We show that the dark tracks are plasma vacuums generated by the 
bouncing and interfering of shocks and expansion waves, upstream an initial slow magnetoacoustic 
shock produced by a localized deposition of energy modeled with a pressure perturbation. The same 
pressure perturbation produces a transverse to the field or perpendicular magnetic shock giving rise 
to nonlinear waves that compose the kink-like plasma void structures, with the same functional 
sunward decreasing phase speed and constancy with height of the period, as those determined by 
the observations. 
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I. INTRODUCTION 



Dark sunward sinuous lanes moving along a fan of rays 
above post-flare loops towards a supra-arcade with a de- 
celerating speed in the range of ~ (50 — 1000) k m s" 1 
were firstly described by McKenzie and Hudson jl999l |. 
using the Yohkoh Soft X-ray Telescope. Further observa- 
tions showed that the fraction of solar flare long-duration 
event arcades possessing fans is not small, and that down- 
flows above such arcades are not uncommon (McKenzie 
2000). The phenomena wer e exte nsively studie d by s ev- 
eral authors ( Wang et al. Il999l ; Innes et al. l2003k .b: 
Svestka et a l. Il998t Asai et al. 12004 McKenzie and 
Savage [20091 ) and different theoretical explanations have 
been offered to give acco unt of the physics of the under- 
lying process. McKenzie |2000j | dismissed the interpreta- 
tion that X-ray voids could be explained as evidence of 
blobs of cool material separated from the main body of 
a CME falling gravitationally back to the surface of the 
Sun. This author supported the explanation that flux 
tubes linking an above current sheet retract into the ar- 
cade u nder the force of magnetic tension (see also Wang 
et al. 1999). The down moving structures observed at 
40 — 60 Mm heights above the top of the arcade indi- 
cated that the reconnection site must be lo cated above 
the flare arcade. Sheeley and Wang 2002j | interpreted 
the phenomenon as the top of collapsing loops with void 



tails in the wakes where reconnection occurs. Asai et 



al. [2004j proposed, based on RHESSI observations, that 
dark lanes are sunward voided flows generated by the re- 
connection process developed by a current sheet at the 
supra-arcade. 

Verwichte et al. ( 20051 hereafter VNC) analyzed trans- 
verse to the magnetic field oscillations associated with 
sunward dark lanes in the post-flare supra-arcade of 
NOAA active region 9906 on the 21st of April 2002. 
They interpreted two oscillating in phase edges (C and 
D in VNC) as fast surface magnetoacoustic kink wave 
trains guided by the ray structure. They founded that the 
phase speeds and the displacement amplitudes decrease 
as they propagate downwards and derived a function to 
fit the amplitude of the oscillating perturbation at a given 
height. The period remained constant with height. As a 
function of time, the phase speeds decreased and periods 
increased. 

In this Letter, based on the following arguments, we 
propose a new scenario to numerically give account of 
the phenomenology described in VNC: 1) the observa- 
tional phase speed magnitude in VNC is several hun- 
dreds of km s _1 lower than the Alfven estimations, this 
could suggest that slow speed sunward processes, rather 
than fast ones, are at work; also, 2) an impulsive en- 
ergy deposition or pressure pulse can lead to lower speed 
void features obtained as a result of nonlinear interact- 
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ing slow shocks, i.e., where the magnetic field plays the 
role of being the wave guide of a funda mentally hydro- 
dynamic shock (Fernandez et al. l2009h . This fact is a 
starting supposition to simulate the dynamics along the 
magnetic field lines. 3) The sunward dynamics is differ- 
entiated, from the transverse one because conductive en- 
ergy transport along field lines is far more effective than 
across them, thus, we assume independence between the 
dynamics of two main directions: the sunward and the 
transverse to the magnetic field one. 4) As the impul- 
sive triggering mechanism acts in all directions, a per- 
pendicular magnetic shock can develop in the transverse 
direction. In conclusion, our modeling is accomplished 
supposing that a pressure pulse applied to the coronal 
inhomogeneous medium, structured by a main magnetic 
field, is a triggering mechanism of two fundamental in- 
dependent shock schemes, an hydrodynamic one and a 
magnetic perpendicular one. On each direction the in- 
teraction of the nonlinear waves will give rise to the dif- 
ferent patterns we are interested in. We emphasize that 
if the simulations reproduce the observational data, e.g., 
the kink-like features, the temporal cadences, the phase 
speed and sunward pattern functional decreases, in the 
frame of a few reasonable suppositions, this implies that 
the overall scheme could be at the bases of the underlying 
physical processes. Even when a full numerical simula- 
tion is required, this two ID ^ modeling, is a useful first 
approach, certainly, if due to its simplification and to the 
accuracy in the observational description suggests that it 
had captured important physical processes. 



II. NUMERICAL METHOD 

To simulate time dependent MHD flo ws we use a 
method developed by us (Elaskar et al. l200lh . The 
equations used are the continuity, momentum, energy 
and magnetic induction equations together with the state 
equation forming a system of partial hyperbolic differen- 
tial equations. The numerical approach consists of an 
approximate Riemann solver with the TVD sche me pro - 
posed by Yee for gasdynamics flows (Yee et al. Il985l) . 
The "eigensystem" technique presented by Powell [l995j | 
was also used with the ei genvec tor normalization pro- 



posed by Zarachay et al. 1994]. The accuracy of the 



technique was verified by means of simulations of the 
Riemann problem introduced by Brio and Wu and the 
Hartmann flows (Maglione et al. l2007t ). 



III. RESULTS AND DISCUSION 

In VNC, to study supra-arcade ray structures the au- 
thors chose a rectangular subfield of parallel (x) and per- 
pendicular (y) to the supra-arcade ray coordinates where 
the structures were visible. They could describe the os- 
cillatory motions of four dark transverse displacements 
for a fixed vertical coordinate as a function of time, i.e. 



x — x a = 45.6Mm. Also, they displayed the inten- 
sity variations of the ray parallel coordinate as a function 
of time i.e., I x (t) averaged over a segment of the y co- 
ordinate, Ay ~ 4Mm, (Fig. 1-2 in VNC). They found 
that the period r remains constant with height but the 
amplitude and phase speed decrease as the wave trav- 
els sunwards. Assuming a transverse perturbation of the 
form 6£y(x,t) = ^4cos[^ — 0], t is the time, and x the 
height, the functions that fit the amplitude and the phase 
speed are respectively 

A{x) = A exp[^i] V ph = V pho (l - 

La V L v 

with correspondent parameters given in Table 1. 

To reproduce the observations we schematically chose 
two set of characteristic values in VNC: edges C and D, 
which seem to have different structure. Edge C's tail has 
a kink-like structure and edge D, more irregular, resem- 
bles a sausage-like structure in part of its tail. Table 1 
shows the observational wave parameters in VNC used 
to adjust the numerical simulation. 

To perform the numerical simulation the physical 
quantities are expressed in a dimensionless form. As 
in VNC the coordinate x represents the sunward direc- 
tion and the y coordinate the transverse to the magnetic 
field one, as seen by the line of sight. Figure [1] shows a 
scheme of the numerical domain of the two 1D^ indepen- 
dent integrations performed, i.e. two main initial shock 
wave fronts triggered by the same pressure perturbation. 
One directed sunwards, parallel to the field lines, and the 
other directed transversally to them. The composition of 
the two processes will give the overall description. The 
details of the energy deposition mechanism are not taken 
into account here e.g., re conne ction processes as the one 
suggested by Asai et al. [2004 1 . 



A. Transverse perpendicular magnetic shocks 

Traditionally it is expected that the observations are 
due to the direct action of boundary effects between the 
void and the surrounding medium. Our simulations show 
that the observational void patterns can also be obtained 
as a result of nonlinear evolution of waves triggered by 
an instantaneous pressure pulse. It is required, however, 
that the media has density inhomogeneities aligned with 
the magnetic field, in order to produce the interferences 
upon reflection that are responsible of sustaining the void 
lanes. The initial perturbation, seen as a transverse to 
the field line phenomenon, produces opposite magnetic 
shock waves that bounce in the lateral and denser bound- 
ary medium composing an oscillatory pattern of nonlin- 
ear wave interactions whose evolution leads to an internal 
defined central vacuum region. The effect of the density 
inhomogeneities -which produce partial reflections and 
absorptions of waves, e.g, it was estimated that density 
values can rise at least between two and five times from 
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the v oid to the denser medium (McKenzie and Hudson 
1999; Inncs ct al. 12003ft - are taken into account as density 
steps modeled by reflecting boundaries at the extremes of 
the integration domain. We performed several simulation 
with different values of the non-dimensional boundary 
density, allowing partial or total rebounds of the waves, 
and we found that the result is robust in the sense that 
the oscillating pattern is conserved for a density range 
going from four to ten times the non-perturbed initial 
value. A study of sensitive to the initial and boundary 
conditions will be matter of a new work together with 
a ID and 2D simulation. A resulting simulation that 
accurately fits the edge C description in VNC is shown 
in Fig. [2^,. The non-dimensional boundary density is 
assumed constant and equal to 5. The color palette indi- 
cate differences in the density values, a partial rebound of 
half of the momentum occurs at y n = and y n = 1. The 
darker region corresponds to the resulting vacuum zone 
pattern, with vanishing values of the density, i.e. the 
non-dimensional vacuum density is 0.3 while the exter- 
nal medium value results ~ 1.4. Note that the simulation 
assumes an initial pressure pulse that triggers the whole 
pattern, thus the following process of wave interactions 
produces and sustains the void cavity for a time interval 
which is in accordance with the observations, as it can 
be seen from the figure. 

If the yi partition is modified, a sausage-like descrip- 
tion can be obtained as it is shown in Fig.^s. This figure 
was adjusted with the sausage-like edge D data. The 
differences between Fig. [2^, and Fig. [2p are due to the 
fact that in the first case an asymmetric rebound occurs 
while in the second one the rebound is symmetric. Kink- 
like structures seem more probable than the sausage-like 
ones because the symmetric condition requirement of the 
external inhomogeneous media is more improbable. Fig- 
ure [5] shows that the vacuum is immediately originated 
by the pulse and the interaction of the non-linear waves 
composes and sustains a kink-like or a sausage-like pat- 
tern for a large period of time. 

The procedure of convertion into the dimensionless 
quantities and the relations that couple the transverse 
and radial parameters are described in next subsection. 



B. Calibration 

The non-dimensional quantities are the ratios of a di- 
mensional quantity and a reference one. We assume a 
reference value of the magnetic field and the gas pressure 
of B re f — 10 G and P re f — 1.38 Pa respectively. To 
derive the dimensional period, amplitude, density and 
speed we proceeded as follows. I) Using the first two 
lines of non-dimensional values in Table 2, we obtained 
Fig. O II) From the figures we determined the mean 
non-dimensional speed chosen as: 
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FIG. 1: Scheme of the loop description. 

where r„ and A n are the non-dimensional periods and 
amplitudes respectively. The correspondent dimensional 
speed is: 

OA 2A 
V = V n V ref = ► V ref = — — (3) 

T T V n 

where the dimensional values r and A are taken from 
Table 1 as reported in VNC for C and D edges. From 
cq. ^ and cq. ([3]), taking into account that V re f — 
Bref /y/l^oPref, we obtain the reference value for the den- 
sity p re f = 1.67 10" 11 Kgm~ 3 . Note that as r = r„r re / 
and that from Fig. [2^,, r„ = 4.5 (or equivalently r n = 2.2 
from Fig. using the edge C value r = 134 s, taken 
from Table 1 (or equivalently r = 175 s from edge D) it 
results: r re / = 29.7 s for the kink structure and r re f = 
79.5 s for the sausage one. Using the same argument for 
the amplitude, A — A Q — A n A re f, with A Q taken from 
Table 1 and A n from Fig. [5] we obtain A re f = 6040 km, 
for the kink structure, and A re f = 2460 km for the 
sausage one. These quantities are also used as the refer- 
ence values for the simulation in the sunward direction, 
together with the initial conditions shown in the last line 
of Table 2. 



C. Radial hydrodynamic shocks 

In the radial direction the same pressure pulse gener- 
ates shocks traveling sunwards and outwards, along the 
magnetic field lines. The nonlinear interaction of out- 
ward rebound waves and sunward absorbed ones com- 
pose an overall descending void structure. The pressure 
pulse is supposed to be associated with a reconnection 
event implying that field lines should retract away from 
the reconnectio n site under the force of magnetic ten- 
sion (McKenzie l2000h . Thus, we suppose that a par- 
tial rebound occurs upwardly near the reconnection site. 
Downwardly, the background initial density is assumed 
constant and equal to the reference value. The sunward 
perturbation is absorbed by the non-perturbed media as 
it travels, however, a total absorption constraint is as- 
sumed at the sunward boundary. Our origin of coordi- 
nates in the radial direction is located at x a = 45.6Mm. 
In VNC x a is the radial coordinate indicating the fixed 
position from where the motion is described, and Table 
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FIG. 2: a) Non-dimensional kink-like structure as edge C in VNC. Asymmetric rebound of the lateral shock in the denser 
media. Dimensional values obtained multiplying by 6040 km and 29.7 s. b) Non-dimensional sausage-like structure as edge D 
in VNC. Symmetric rebound of the lateral shock in the denser media. Dimensional values obtained multiplying by 2460 km 
and 79.5 s. 



1 values were measured. We chose Xt — 50Mm as the 
position of the triggering pressure pulse. Thus, Xt is also 
the distance traveled by the void features, i.e., taking 
into account different trajectories of the void heads in 
VNC observations this distance results a good assump- 
tion. The non-dimensional value of Xt is x on = 0.7 and 
the non-dimensional initial density, p n — 1, corresponds 
to the dimensional p re f value. Figure [3^ is equivalent 
to Fig. (2C) in VNC. Note that the radial perturbation 
initially generates two shocks, one moving sunwards and 
the other upwards. Figure [3J3 is a zoom of Figure [3^, 
showing the upward moving features, which are forced to 
partially rebound at x n = 1. The upwardly moving fea- 
tures in Fig. [3Jd correspond to a shock front, recognized 
because of the density enhancement at the shock front, 
and a contact discontinuity wave, recognized because the 
pressure and the speed of the flow are not changed while 
the wave passes. Also, an expansion wave travels down- 
wards, it is recognized because the density is diminished 
at its pass. Figured shows that the interaction of sub- 
sequent waves composes a downward moving void as it 
is indicated in by the color contrast, i.e. darker features 
correspond to lower values of the density. The clearest 
curve that originates at x n — x on — 0.7, t n = 0, and 
ends at x n — 0, t n — 0.9, is the main downward shock, 
i.e. while crossing the shock front the density varies from 
p n = 1 to p n = 3. Later interactions compose a region 
of smooth expansion that moves sunwards until it is ab- 
sorbed at x n = 0, t n ~ (0.9 — 2.15). Further interactions 
generate a low density upper zone, also formed by the 
coupling of nonlinear upward and downward moving per- 
turbations, it is separated from the later zone by a sharp 
curve. The variation of the density as a function of the 
sunward direction for fixed times can also be seen from 



Edge 


t(s) Vphoikms^ 1 ) 


L v (Mm) 


A (km) 


L A (Mm) 


C 


134 202 


40.4 


906 


34 


D 


175 176 


14.3 


246 


11.5 



TABLE I: Observational wave parameters taken from VNC. 



Fig-Et- For the non-dimensional time t n = 0.594 we note 
a first peak corresponding to the main shock front, the 
posterior smooth decaying values of density corresponds 
to the expansion region. The interaction with descend- 
ing perturbations determine another shock peak with a 
descending step of density at ~ 0.41. The behaviour is 
similar for t n — 0.95 without the second shock. At larger 
times, t n = 1.38, as the main shock was already absorbed, 
only the step of the density, while crossing the voided re- 
gion, is registered. 

To compare our speed values with those in VNC we cal- 
culated the speed of the curve limiting the vacuum zone. 
The VNC curve was obtained with the speed function of 
cq. ([T]) replacing x by 0, in accordance with our choice 
of the origin of coordinates. Figure 0] shows the com- 
parison between the numerical and observational curves. 
Note that the the reference value V p ho — 202kms _1 , the 
speed measured by VNC at x , is the simulation speed 
at the origin of coordinates. The figure shows the good 
correspondence between the observational velocity curve 
and the numerical one. 



IV. SUMMARY 

We have integrated the MHD ideal equations to simu- 
late observational dark lane data obtained by Verwichte 
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FIG. 4: Numerical expansion curve vs. observational VNC 
data. At x = x a , origin of coordinates, V p h = V p ho 
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FIG. 3: (a) Non-dimensional sunward displacement as a func- 
tion of x n and t„ for edge C. Dimensional values obtained 
multiplying by 48.6 Mm and 29.7 s. Equivalent to Fig. (2C) 
in VNC. (b) Zoom of (a) showing details of the upper re- 
bound, (c) Density profile as a function of height for the 
times t = 0.594, 0.95, 1.38. 



Edge 


yi,V2,y3 


P 2 /Pl,3 


Xl, x 2 ,x 3 


C 


0.7,0.04,0.26 


10 




D 


0.48,0.04,0.48 


10 




C(radial case) 




10 


0.7,0.25,0.05 



TABLE II: Initial condition values for the transverse simu- 
lation shown in Fig. [2] first two lines. Initial condition val- 
ues for the sunward simulation shown in Fig. [3] third line. 
Pi = Ps, pi = P2 = pz = 1. 



et al. [2005} (VNC). We simulated the effects of an initial 
impulsive and localized deposition of energy -supposed 
to be associated with above reconnection processes- in 
a plasma structured by sunward magnetic field lines. 
The impulsive phase is modeled by a pressure perturba- 
tion that initiates two main different type of processes, 
a fundamentally hydrodynamic shock pattern directed 
sunwards and a perpendicular magnetic shock one, i.e., 
transversal to the magnetic field. The two patterns are 
supposed to be independent processes however linked by 
their common origin and background magnetic and den- 
sity conditions. Firstly, we reproduced the oscillatory 
kink-like pattern of the transverse y component. Vary- 
ing the initial conditions we were able to calibrate and 
reproduce the C edge in VNC. Also a sausage-like pat- 
tern was obtained resembling the tail edge D in VNC. 
The differences between patterns are due to the non- 
symmetric (most probable kink-like) or symmetric (less 
probable sausage-like) rebounds of the magnetic trans- 
verse shocks in the denser external medium. The result- 
ing interactions of bounced nonlinear shocks composes 
and sustains, in accordance with observational character- 
istic times of the phenomenon, a density structure with 
a central void resembling a kink-like mode or a sausage- 
like mode as shown in Fig. rj^-b. 

The same reference values were used as initial condi- 
tion to simulate the sunward evolution. We could qual- 
itatively reproduce the observational data showing that 
initially two opposite shock evolving fronts are produced. 
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One evolves towards the sun surface until is absorbed, 
and the other is forced to rebound upwardly resembling 
the action of a reconnection site. Afterwards, the in- 
teraction of upwardly and downwardly moving pertur- 
bations form an expansion wave region that lowers the 
density of the medium. Also a voided zone is formed 
and sustained, lasting for times comparable with obser- 
vational times, due to the continued interfering of non- 
linear waves. This evolution is described in Fig. [3] which 
reproduces with good accuracy Fig. (2C) in VNC. Fi- 
nally, the composition of the two dynamics, one corre- 
sponding to the sunward directed hydrodynamic shock 



pattern and the perpendicular magnetic shock, composes 
an overall description of a moving transversally shaking 
void, resembling a kink-like mode as it moves towards 
the sun surface. 
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